clear all
set scheme s2color
cd ..
use "Data/Processed Data/ACSdata_env.dta" 
lpoly env_ben_PC Median_HH if Median_HH<150000, kernel(ep) bw(10000) pwidth(15000) noscatter ci title("") xtitle("Median Household Income") ytitle("Environmental Benefits ($) per Capita") legend(label(1 "95% CI") label(2 "Fitted Env Benefits") position(11) ring(0) col(1)) note("") graphregion(color(white)) bgcolor(white)
graph export "figures\env_ben_vs_income_wsize.pdf", replace
